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Abstract: The spectral energy distribution and variability of several blazars (Mkn 501, 1 ES 2344+51.4, PKS 2155-30.4, 
1 ES 1218+30.4, 3C 454.3) along the blazar sequence is investigated using a selfconsistent and timedependent lepto- 
hadronic hybrid emission model. The BL Lac objects in the list are successfully modelled with purely leptonic radiation 
processes (essentially Synchrotron Self-Compton emission), while the Flat Spectrum Radio Quasar requires highly rel- 
ativistic hadrons to be present within the jet. Variability is exploited as well with our model to distinguish between 
Self-Compton and hadronic gamma radiation making use of their different signatures in lightcurves. As a consequence 
active galactic nuclei with 2; > 0.5 are the best candidates as sources of extragalactic consmic rays, since High-Peaked 
BL Lac objects do not seem to accelerate protons to energies above thermal. Furthermore the parameters found during the 
modelling of the objects agree very well with obervations of e.g. superluminal motion or typical variability timescales. 
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1 Introduction 

Blazars have gained a lot of attraction in the past decades, 
mainly due to the discovery of their very high energy 
(VHE) emission up to 10^^ Hz with Air-Cherenkov Tele- 
scopes like H.E.S.S. or MAGIC. Ever since more and more 
blazars with a redshift up to z = 0.536 (3C 279 with 
MAGIC |6|) are discovered in the VHE making a system- 
atic investigation possible. The spectral energy distribution 
of blazars shows two prominent humps, one in the optical 
to hard X-Ray band and one in gamma-ray energies. The 
emission is variable on short time scales across the whole 
spectrum with large amplitude changes. Fermi extended 
our view allowing a nearly all-sky survey between 20 GeV 
and 300 GeV. Thus the availability of multiwavelength 
(MWL) data, even during outbursts of blazars, is rapidly 
increasing. 

Blazars are phenomenologically divided into different sub- 
classes based on their peak frequency forming the blazar 
sequence: Flat Spectrum Radio Quasars (FSRQs) like 
3C 279 as the most luminous objects having their first peak 
in the near infrared to optical regime. As the overall lumi- 
nosity decreases in Low frequency and Intermediate fre- 
quency peaked BL Lac objects (LBLs and IBLs respec- 
tively) the first hump shifts to higher energies with high 
frequency peaked BL Lac objects (HBLs) achieving the 
highest photon energies up to 2 TeV showing the lowest 
luminosity among blazars ifTTl 1241 . 

While HBLs are well explained by Synchrotron Self- 
Compton (SSC) models [[T6ll26l . i.e. the first peak is due to 



synchrotron radiation of primary electrons and the second 
hump arises from Compton upscattering (IC) of the very 
same radiation by the primary particles, LBLs and FSRQs 
may not be explained the same model. Their hard VHE 
radiation and wide spread of the peaks are not consistent 
with IC scattering of synchrotron photons IfTOl . The ori- 
gin of the VHE emission in this subsample is still under 
passionate debate: either external radiation providing the 
target photons for IC upscattering within the jet lfT4l . or it 
is due to photo-hadronic cascades initiated by high energy 
protons combined with their synchrotron radiation lfT0l[T9ll . 
In this paper we use a selfconsistent, timedependent hybrid 
model co-accelerating electrons and protons with Fermi-I 
and Fermi-II processes at the very same site while taking 
all the relevant loss processes into account to reproduce the 
emission of different blazars along the sequence with just 
one model by varying the parameters like the p+ to ra- 
tio. Outbursts along the sequence can be exploited with the 
model as well, supporting e.g. the discrimination against 
EC enhanced models or the determination of sources ac- 
celerating protons and not only electrons to the highest en- 
ergies. 

Consequently FSRQs are the most promising candidates as 
sources of extragalactic ultra high energetic cosmic rays 
(UHECRs) among blazars. The UHECRs are believed to 
be extragalactic since we clearly see the GZK-Cutoff in the 
well known spectrum of CRs lU |4l and due to the lack of 
a large-scale anisotropy in the spectrum towards the galac- 
tic disk. Of course these UHECR statements are strongly 
limited by statistics. 
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2 The Model 

In this section we briefly introduce our model focusing on 
the additions compared to the more detailed description al- 
ready available in |25 , 26 1. The geometry and general prop- 
erties of the model stay untouched while the corresponding 
one dimensional Vlasov-equation in the diffusion approx- 
imation (see e.g. 1231 ) is solved for electrons and protons 
separately in two spatially separated zones. The upstream 
spherical acceleration zone is nested in a bigger radiation 
zone. The kinetic equations in the acceleration zone there- 
fore are 

dtn, =d^ [{f3s,a^ - t-^la) ■ ^i] + 

[[{a + 2)U,,:,]-'j^d^n,] + Qo,. - (1) 

''esc,i 

with i being electrons (e) or protons (p) respectively. The 
synchrotron /3s and iesc depend on the considered particles' 
mass as follows (CGS units) 
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By setting iacc/^esc as a parameter during the modeling the 
acceleration timescale will scale linearly with the mass as 
one would expect, i.e. 7max.i particles can reach scales with 
mf, whereas the spectral index being independent since the 
acceleration mechanism itself stays the same, i.e.: 
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Note that the only new parameter introduced compared to 
the simple SSC case is the injection height Qo.p of protons, 
but due to the opacity to p + 7 processes this has a major 
effect on the radiation zone. The plasmaphysical properties 
like i^Tij, Wshock or UAifven can still be calculated as described 
in|l25]|26l. 

Protons with non-thermal energies are above the threshold 
for pion production p + -y p + ngTr'^ + n+7r+ + n_7r^ 
within the jet which decay into electrons, positions and 7s 
(and neutrinos), the channels are 
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This will result in the third contribution, besides proton 
synchrotron and IC from electrons, to the VHE peak in the 
SED of blazars when hybrid models are considered. We 
use the Monte-Carlo simulations of ifTSl to determine the 
production rate of the stable electrons, positrons and 7s 
from the interactions of eq. (|4|. We do not account for 
synchrotron losses of the unstable tt*, i.e. the photo- 
hadronic interactions are assumed to be instantaneous com- 
pared to the synchrotron loss timescale. The produced 7s 
have energies well above 511 keV and will therefore result 
in pair cascades of e*. Other processes are not taken into 



account since they play a minor role in typical AGN jets 
Q. 

The kinetic equation of electrons in the radiation zone be- 
comes 

dtN,^ =d^ [(/3,,e7' + ^^10(7)) -Ne-]- 

Qpp + Qp.- + fc'^ (5) 
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with b = Racc/Rrad and the pair production rate Qpp(7) 
calculated from ifTTI . The production rate of e~ from tt~ 
is calculated using fTSl eq. (30) with the corresponding 
^-function of their paper. For positrons one finds 



dtN,+ =d^ [(/3,,e7' + ^IC(7)) • A^e+] 
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(6) 



since there are no primary positrons in our model. Comp- 
ton scattering losses of protons as well as p + 7 losses can 
be neglected following the discussion of (|9][H], hence 



rad,esc,p 



(7) 



The Qpp, etc. gains and IC losses of course depend on 
the photon field which is calculated from the corresponding 
production and loss rates: 



dtN^ ^ Rs + Rc + R^o ~ c (assA + ctpp) 



(8) 



with the additional app accounting for losses due to the 
pair-channel calculated from lfT3l . eq. (4.3) and the 
vr'^-decay channel as calculated from fTSl, eq. (11). The 
synchrotron production rate Rs now consists of electron, 
positron and proton synchrotron radiation being produced 
due to eqn. (|5]l to Q, Rc takes IC photons from and e+ 
into account. 

The model is tested with respect to energy and particle con- 
servation of all contributing processes. Note that we do not 
introduce any ad hoc assumptions and all spectra arise self- 
consistently during the time evolution of the model. This 
allows conclusions on the microphysics of the jet exploit- 
ing the variability of blazars as well, see [26|. 



3 Results 

The blazars used in the modelling summarized in Ta- 
ble [T] show variability from the sub-hour timescale 
(PKS 2155-30.4) to outbursts lasting for days (Mkn501, 
1 ES 1218H-30.4) or at least have been observed in multi- 
wavelength observation campaigns in different flux states 
(e.g. 1 ES 2344H-51.4). Table [1] also shows the the model 
parameters found for each object for its lowstate emis- 
sion as well as their blazar-type and redshift. All pa- 
rameters displayed in Table [T] are in a physically sensi- 
ble range compared to other observations like superluminal 
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Source 


z 


Type 


Qe(cm"^) 


Qp(cm"^) 


B{G) 


Race 


(cm) 


^rad 


(cm) 


^acc/^esc 


a 


6 


Mkn 501 


0.034 


HBL 


3.1 • 10^ 





0.09 


2.2 • 


10" 


5.0- 


10^^ 


1.20 


50 


34 


1 ES 2344+51.4 


0.044 


HBL 


6.5 • 10^ 





0.10 


1.8 • 


10" 


6.0 • 


1015 


1.05 


50 


18 


PKS 2155-30.4 


0.117 


HBL 


8.0 • 10^ 





1.4 


1.0 • 


1013 


5.0- 


1014 


1.13 


20 


49 


1 ES 1218+30.4 


0.182 


HBL 


6.3 • 10" 





0.12 


6.0 • 


10" 


3.0 • 


1015 


1.11 


10 


44 


3C 454.3 


0.859 


FSRQ 


3.8 • 10^ 


4.2 • 10^ 


10.2 


5.0 • 


10" 


5.0- 


1015 


1.10 


5000 


43 



Table 1: List of blazars used along with the parameters used for modelling the lowstate emissions with our model, see text 
for details. 



motion putting upper limits on S or the observed variabil- 
ity timescales / Schwartzschlid radius lower boundaries on 
i?iad and coincide within an order of magnitude with other 
modelling attempts by e.g. F. Tavecchio (1 ES 2344+51.4, 
see |22|) while not relying on ad-hoc assumptions. The 
magnetic field for each blazar is high enough to con- 
fine the particles within the emission region, i.e. Tgyi « 
jniic/eB ^ i?rad- Without introducing a significant num- 
ber of non-thermal protons Qp in the case of 3C 454.3 one 
fails to find parameters modelling the SEDs satisfying the 
mentioned criteria hence strongly disfavoring a purely SSC 
scenario for FSRQ like objects. 

With our timedependent code we are in the unique posi- 
tion to exploit the variability of blazars from HBLs to FS- 
RQs providing information about the nature of the physics 
within the jet. This can be used to i) constrain the model 
parameters more tightly allowing together with the self- 
consistency of the model microphysical interpretations of 
them and ii) to discriminate between purely leptonic jets 
and hadronic scenarios since each leaves a typical imprint 
in lightcurves. 

For the HBLs 1 ES 1218+30.4 and PKS 2155-30.4 the 
timedependent results can be found in |25| and |26l re- 
spectively. The lightcurves in the VHEs can be traced re- 
markably well simply by injecting more primary electrons 
into the blob. In the case of Mkn 501, as well being mod- 
elled purely leptonic, the variability seems to have a dif- 
ferent nature which can clearly be seen in Fig. [T] Unlike 
1 ES 1218+30.4 or PKS 2155-30.4 the peak in the VHEs 
shifts towards higher energies during the outburst observed 
by VERITAS and Fermi in May 2009. Similar behaviour 
was, even more prominently, observed by BeppoSAX dur- 
ing the famous outburst 1997 in the X-Rays ||2T1 making 
the scenario somehow persistent even when considering 
two unrelated events (blobs). The recent flare was modelled 
as a linear dropoff in the magnetic field from B = 0.09 G to 
B — 0.043 G over a time of tdrop = 11 d which consistently 
is accompanied by a dropoff in the injection of particles to 
Qe = 1.60 ■ 105 cm~3 and of iacc/^rad to 1.0. After another 
^high = 40 d these values linearly took their lowstate values 
displayed in Tab. [T]within irise — 28 d. The average over 
the whole outburst results in the highstate (dashed line in 
Fig. [TJ of Mkn 501, all times are given in the observer's 
frame. The highstate of 1 ES 2344+51.4 was found by set- 
ting (5 = 23 and Qe = 3.75 • 105 cm^'^ instead of the values 
found in Tab. [T] suggesting observing radiation from two 
different blobs, therefore no time evolution was tested at 
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Figure 1: Modelled lowstate emission of Mkn 501 (solid 
line, parameters see Table [U and the averaged outburst in 
05/2009 (dashed line, averaging and parameters see text). 
The lowstate was found taking the multiwavelength obser- 
vations of |3| as well as |5| (brown points) into account. 
The VHE data has been EBL corrected, see |i3J. 

this particular HBL, although a speed-up or altering angle 
to the line of sight of the emitting region might be possible 
but unlikely interpretations. 

Concerning 3C 454.3 a significant number of protons had 
to be injected into the blob which are accelerated due to 
Fermi I and Fermi II processes up to 10^ GeV for a small 
fraction of them. Purely leptonic fits lead to arbitraryly 
low magnetic fields as well as S ^ 100 and thus were 
neglected. The radiation in the optical regime essentially 
is synchrotron emission from primary electrons, in the X- 
Rays the SED is dominated by proton synchrotron and re- 
processed radiation (synchrotron radiation of pairs cas- 
cading down from pion production and their radiation), 
see Fig. [3] In the Fermi-LAT energy range the SED of 
3C 454.3 seems to be dominated by proton synchrotron 
emission, see Fig. [3] Concerning variability the proton 
synchrotron radiation in the X-rays is providing a very sta- 
ble flux-basis with the proton synchrotron timescale of the 
order of months. Shorter timescales as observed during the 
flares in December 2009 |20| are hence due to cascade ra- 
diation being enhanced by fluctuations in the primary elec- 
trons and/or their synchrotron emission itself due to B field 
changes (the latter might explain the observed flaring be- 
haviour of I.20J ). In both cases the cascade radiation is out- 
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Figure 2: Modelled SED of 1 ES 2344+51.4 in two differ- 
ent flux states (lowstate: solid line, parameters see Table [T] 
highstate: dashed line, parameters see text). The data was 
was almost simultaneously taken, see I.22J for details. 

matching the proton synchrotron radiation of the lowstate 
while cascading down. 

4 Conclusion and Discussion 

We were able to model all four low-redshift blazars con- 
sidered with our selfconsistent model without introducing 
non-thermal protons within the jet. Making use of the 
timedependence of our model we constrained the param- 
eters down to the values found in Table [T] However the 
variability of HBLs seems to have no common cause. In 
the case of PKS 2155-30.4 and 1 ES 1218+30.4 inject- 
ing more particles into the blob (i.e. density fluctuations 
within the jet, as the blob travels through them) reproduces 
the high states and lightcurves. Concerning Mkn 501 the 
variability seems to have a different nature: a dropoff in 
the magnetic field accompanied with a consistent dropoff 
in the particle density. This sort of outburst was observed 
in the very recent flare we modelled and also back in 1997 
by BeppoSAX in the x-rays. This might be related to the 
stability of the jet. The variability of 1 ES 2344+51.4 sug- 
gests observing radiation from two different blobs, see 
Sect. 3 for details. Nevertheless all four HBLs allow 
for large amplitute, short time variability and can be ob- 
served at very different flux states over the years as e.g. 
Mkn 501. The situation is different for the rather distant 
FSRQ 3C 454.3. The SED in the Fermi energy range seems 
to be dominated by proton synchrotron radiation which has 
a variability timescale of the order of days as observed 
in outbursts ll20l [TSl . Short time X-ray variability is ex- 
plained by cascade radiation dominating over the promi- 
nent lowstate proton synchrotron emission during an out- 
burst. The parameters found for 3C 454.3 agree well with 
other hadronic models of by M. Bottcher fTOI for other FS- 
RQs, e.g. 3C 279. 

Concluding HBLs do not allow hadrons to be accelerated 
to energies above thermal due to their relatively low mag- 



Figure 3: Modelled SED of 3C 454.3 including flieir indi- 
vidual processes, dashed red: proton synchrotron, dashed 
blue: subsequent cascade radiation. The grey line shows 
the modelling attempt using only er with physically mo- 
tivated parameters which fails to reproduce the VHE part. 
The data was taken from |2|, the model parameters can be 
found in Table [T] The VLT observations are not simultane- 
ous. 

netic field being not able to confine them within the emit- 
ting region for a significant time. Hence FSRQs seem to be 
the most promising candidates among blazars as sources 
of extragalactic cosmic rays. However, their lower peak 
energies compared to HBLs point towards maximum pro- 
ton energies (after accounting for the blob's bulk motion) 
above the ankle in the spectrum of CRs but not to ener- 
gies of UHECRs. Acknowledgements MW thanks Elitenet- 
zwerk Bayern and GKl 147 for their support. FS receives 
support from DFG through grant SP 1 124/1. 
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